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We numerically investigate the radial dependence of the velocity and temperature fluctuations 
and of the time-averaged heat flux j(r) in a cylindrical Rayleigh-Benard cell with aspect ratio F = 1 
for Rayleigh numbers Ra between 2 x 10^ and 2 x 10^ at a fixed Prandtl number Pr = 5.2. The 
numerical results reveal that the heat flux close to the side wall is larger than in the center and that, 
just as the global heat transport, it has an effective power law dependence on the Rayleigh number, 
j(r) oc Ra'^'^'^^K The scaling exponent 7j(r) decreases monotonically from 0.43 near the axis (r ~ 0) 
to 0.29 close to the side walls (r ~ D/2). The effective exponents near the axis and the side wall 
agree well with the measurements of Shang et al. (Phys. Rev. Lett. 100, 244503, 2008) and the 
predictions of Grossmann and Lohse (Phys. Fluids 16, 1070, 2004). Extrapolating our results to 
large Rayleigh number would imply a crossover at Ra ~ 10^^, where the heat flux near the axis 
would begin to dominate. In addition, we flnd that the local heat flux is more than twice as high 
at the location where warm or cold plumes go up or down, than in the plume depleted regions. 

PACS numbers: 47.55.P-,47.55.pb,44.25.+f 



I. INTRODUCTION 

In Rayleigh-Benard (RB) convection a fluid is heated 
from below and cooled from above. This problem of ther- 
mal convection is of the utmost importance from an ap- 
phed point of view. Examples are thermal convection 
in the atmosphere, in the oceans and in process tech- 
nology. For recent reviews of RB convection we refer to 
Refs. [HH]. 

In a cylindrical container, the dynamics of a RB system 
depends on three control parameters, the Rayleigh num- 
ber Ra = g/3AL^ /i^K, the Prandtl number Pr = v/k, and 
the aspect ratio T = D/L. Here g is the gravitational ac- 
celeration, /3 the isobaric thermal expansion coefficient, 
A the temperature difference between the top and bot- 
tom plates, V the kinematic viscosity, n the thermal dif- 
fusivity, and L and D are the height and diameter of the 
cylinder. The response of the system is expressed by the 
Nusselt number Nu, the dimensionless heat flux [1]. 

Previous studies mainly focused on determining the 
global heat flux as a function of Ra and Pr. For water 
(Pr = 4 — 6), in the experimentally available range of 
Ra = 10^ — 10^^, one finds that the global heat trans- 
port effectively scales as Nu ~ Ea^ ^s-o-si [HEHg. The 
effective exponents for the global heat flux are well de- 
scribed by the unifying theory of Refs. [U EHU]- That 
theory also made predictions for the scaling exponents of 
the local heat flux in the center of the cell and at the 
side wall |TT]. The reasoning is based on splitting the 
thermal energy dissipation field into its plume and back- 
ground contributions; similarly the kinetic energy dissi- 
pation is decomposed into its boundary layer and bulk 



contributions. By doing this Grossmann and Lohse 
accounted for the various scalings in the Ra — Pr pa- 
rameter space. They found that the local heat flux has 
an effective power law dependence on the Ra number, 
j(r) cx Ra^^'^^\ and obtained a prediction for the scaling 
exponent 7^ = 0.45 in the bulk (center) and 7^ = 0.22 at 
the plume (side wall) regions. 

In order to understand the heat flux one has to ei- 
ther rely on Eulerian [12] or Lagrangian [131 [13 mea- 
surements where the complex interplay between velocity 
and temperature can be studied. Advancements in the 
experimental techniques made it possible to measure the 
vertical local velocity ^^(r,*) and the local temperature 
r(r, t) at a given spatial location r as functions of time. 
This allowed Shang et al. [TH [TS] to determine the local 
convective heat flux [16] 

^ „.(M)[T(M)-r„| ^ 

where Tq is the mean bulk temperature. They deter- 
mined the probability density functions (PDF) of the lo- 
cal heat flux in the axis and side wall regions and showed 
that the vertical heat flux is highly non-Gaussian and 
intermittent due to thermal plumes. This work stimu- 
lated Ching et al. [T7] to theoretically study the problem. 
They decomposed the velocity field into a part associated 
with strong temperature fluctuations plus a background 
and found that, with the definitions they used, the local 
heat transport associated with the former velocity com- 
ponent near the axis of the cell scaled as i?a^/^. Later 
experiments by Shang et al. [18] revealed that the ef- 
fective scaling exponent for the local heat flux is about 
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FIG. 1. Computational grid used for Ra = 2 x 10^. We 
position numerical probes uniformly in azimuthal direction in 
circles at seven radial locations r/L = 0.06, 0.12, 0.19, 0.24, 
0.34, 0.40, 0.45 (shown in black). For better statistics we use 
two more sets of probes (blue, red) with an offset of one grid 
point. Inset shows a close up view of the probe locations on 
the grid. 



0.49 near the cell axis and about 0.24 near the side wall, 
which confirms the main results obtained by Grossmann 
and Lohse [11], but is in disagreement with the model of 
Ching et al. [T7]. 

In experiments it is quite difficult to measure the heat 
flux at each spatial point in the cell due to problems 
with measurement techniques and the presence of in- 
herent noise levels. The measurements at just one or 
two points may not be enough to understand the com- 
plex dynamics involved in the convection process. The 
present paper offers numerical results which complement 
the work initiated by Shang et al. [TH HHj . We provide 
information on the heat flux at one quarter, one half, and 
three quarters of the cell height for several radial posi- 
tions r, not only near the axis and the side wall. This 
information allows us to understand the two limits in the 
unifying theory |llj on the effective scaling exponents, 
one valid in the bulk (the central region of the cell) and 
the other valid in the plume region. As a result we clearly 
see persistence of the inhomogeneous nature of the flow 
in the radial direction which leads to different scaling 
exponents. 

In experiments the local velocity is measured by 
LDV/PIV techniques at a spatial position which slightly 
differs from the location of the local temperature mea- 
surement. This spatial misalignment may possibly affect 
the results. In numerical simulations, in contrast, one 
has all the information on the flow field and thus the 
local heat flux can be calculated from the velocity and 
temperature measurements at exactly the same position. 



Finally we comment on the local heat flux distribution 
with respect to the large scale circulation and the ulti- 
mate regime mentioned in Refs. [HI HH HH] and provide 
data to illuminate the differences between the measure- 
ments of Shang et al. (18J and the earlier Ching et al. [T7] 
claims. 

In simulations the Ra number range and the duration 
available for time-averaging are more limited than in ex- 
periments. In order to mitigate the latter shortcoming 
we have limited the Ra number of our simulations. In 
the next section we briefly describe the numerical proce- 
dure before discussing our results on the local heat flux 
and the velocity and temperature fluctuations. 



II. NUMERICAL METHOD 

We performed direct numerical simulations for a 
Boussincsq fluid in a unit aspect ratio (F = 1) cylin- 
der with constant temperatures applied at the top and 
bottom plates and an adiabatic side wall. The fluid simu- 
lated in our calculations is water at 32 °C (Pr = 5.2) for 
2 X 10® < Ra < 2 X 10^. The governing equations for mo- 
mentum, energy and mass conservation in dimensionless 
form are given by |20| 
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= -Vp + 9z 
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Pr 
Ra 



1/2 



V^u, 



V^9, V-u = 0. 



(2) 
(3) 



Here the dimensionless variables are the velocity u, tem- 
perature 9 and pressure p (minus the hydrostatic contri- 
bution). The material derivative is denoted by D/Dt. 
The unit vector z is in the direction opposite to gravity. 
The physical variables such as length and velocity are 
made non-dimensional by the cylinder height (L), and 
the free-fall velocity U — y/g/3AL. Constant dimension- 
less temperatures of 1 and are applied at the bottom 
and top plate, respectively. 

The governing equations are solved on a staggered grid 
with second-order accuracy in space and time. For the 
time advancement a third order Runge-Kutta scheme is 
used. This method is stable for a CFL number up to -s/S, 
and here we have limited it to 1.2 by adjusting the time 
step according to the maximum velocity before execut- 
ing each time step . More details about the numerical 
method can be found in Refs. [20l [22j [23] . For the spa- 
tial resolution we followed the criteria set by Stevens et 
al. [24,. 

A summary of the simulation parameters is shown in 
Table H] The flrst and second columns are Ra and the 
number of grid nodes. The volume- and time-averaged 
global heat transport, i.e. the Nusselt number Nu = 

1 + VRaPr 



{u,9), 



is shown in the third column; here 



the overline denotes time averages and the angular brack- 
ets {■)v volume averages. The number of grid points used 



3 



TABLE I. Summary of simulation parameters: The number of grid points used in angular (N^), radial (A^r), and axial 
directions (N^), volume and time averaged Nusselt number (Nu), the number of grid points in the thermal boundary layer 
(nti), convergence of exact relations for ek and et, comparing maximum of grid spacing in angular (Scprn) and axial (Szm) 
directions with the Kolmogorov length scale based on global kinetic energy dissipation rate (77) and the averaging time considered 
for the simulations are shown. The reported times are always measured in the units of free fall time and the lengths in terms 
of cylinder height. For statistical averages we discarded the initial 140 dimensionless times, to prevent transient effects in the 
results. 
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2 X lO** 


193 X 49 X 129 


10.93 


19 


1.008 


0.968 


2.26, 1.45, 3.41 


4200 


1 X 10^ 


257 X 65 X 193 


16.58 


18 


1.010 


0.971 


2.44, 0.73, 2.04 
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2 X 10^ 


257 X 65 X 193 


20.71 


18 


1.007 


0.843 


2.44, 0.73, 1.62 


3700 


6 X 10^ 


321 X 97 X 239 


28.48 


14 


0.989 


0.923 


1.96, 0.97, 1.13 
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1 X 10* 


385 X 129 X 257 


33.25 


15 


0.997 


0.955 


1.62, 0.73, 0.95 
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2 X 10* 
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13 


1.002 


0.937 


1.62, 0.73, 0.76 


2700 
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513 X 161 X 321 
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13 


1.005 


0.947 
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FIG. 2. PDFs of the local instantaneous heat flux j near the 
side wall at r/L — 0.45 (solid, thick lines) and near the axis 
r/L = 0.06 (dashed, thin lines) at mid-height for Ra = 2 x lO'^ 
(red), 1 X 10* (green), 2 x 10^ (blue). The inset shows time 
series for the heat flux near the side wall (thin, red) and in 
the center (thick, gray) for Ra = 1 x 10*. 



FIG. 3. x PDF{j) vs. j is shown for different Ra near 
the side wall region showing convergence of the flatness at 
that location. Note that the y-axis is given in log-scale, in 
contrast to the figures shown in the convergence test by Belin 
et al. [5^, where x PDF{y) is given on a linear scale. The 
inset shows that the flatness F4, increases as a function of Ra. 



to resolve the thermal boundary layers is shown in the 
fourth column. In a RB cell with no slip velocity con- 
dition on the walls the dimensional thermal energy dis- 
sipation rate, St ~ KA^L^^(|V0p)v and kinetic energy 
dissipation rate £k = v^L~'^Pr~^Ra{\Vu\'^)v satisfy ex- 
act relationships with Nu (see Refs. [U |25]), namely 
Et = kA^L'^Nu and Ek = v^L-^Pr-'^Ra{Nu - 1). 
In order to validate our simulations the obtained energy 
dissipation rates are compared with Nu in Table |l] These 
ratios are near one, which proves the adequacy of the grid 
resolution. In the seventh column the largest grid spac- 
ings in the azimuthal and axial directions are compared 
with the Kolmogorov length. In the last column the to- 
tal time used for the statistical averages is shown in free 
fall time units {L/U). The total computational time was 
around 2.2 x 10^ CPU hours on a Power6 computer. 
We placed 1880 "numerical probes" at different ra- 



dial locations on three different horizontal planes (at 
2;/L=0.25, 2:/L=0.50 and z/L=0.75), to obtain point- 
wise data on the temperature and vertical velocity in or- 
der to calculate the local heat flux according to equation 
([1]). In each horizontal plane a number of azimuthally 
nearly equally spaced probes were placed on seven cir- 
cles with radii r/L = 0.06, 0.12, 0.19, 0.24, 0.34, 0.40 
and 0.45, see Figure [T] On each circle we distributed 60 
probes at z/L = 0.50 and 20 probes at z/L = 0.25, 0.75. 
In addition, on the mid-plane, each circle was comple- 
mented by two other circles, one inside and one outside, 
spaced by one radial mesh length as shown in the upper 
right corner of Figure [T] In total, we have information 
from 60 (azimuthal) x 7 (radial) x 3 (sets)= 1260 probes 
at mid-height plane. For the planes at z/L ~ 0.25 and 
at z/L ~ 0.75 we have information from 20 x 7 = 140 
probes. 

As a further check on the consistency of the calcu- 
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FIG. 4. Centered PDFs of the normalized heat flux j* = {j — 
j)ljrma On the mid-plane (a) near the side wall r/L = 0.45, 
and (b) near the axis r/ L = 0.06. 

lations we compare the time- and area-averaged heat 
flux as monitored by all the mid-plane probes mentioned 
before complemented by 340 additional ones uniformly 
spaced in the angular direction and non-uniformly in 
the radial direction. For Ra = 5 x 10® we found a 
global average Nusselt number of 51.58 while from the 
r-averaged heat flux from the probes at mid-height we 
found Nu = 50.83. Similarly, for i?a = 2 x 10^, the two 
results were 78.80 and 76.41. In view of the discretiza- 
tion error introduced by the sparseness of the numerical 
probes (5% of the actual computational cells on the mid- 
plane), this can be considered as a good agreement. 



III. RESULTS 

A. Local heat flux 

In Figure [2] we show PDFs for the local instantaneous 
heat flux j on the mid-plane, at two different radial po- 
sitions, one near the cell axis {r/L — 0.06) and one near 
the side wall {r/L = 0.45). Note that the latter is outside 
the kinematic boundary layer (BL) which at this Prandtl 
number has a thickness of K/L « 3.6i?a""-2^=^° °3 [27]. 
According to this scaling law we have A„/L ss 0.05 for 
Ra = l(f and A„/L w 0.016 for Ra = 10^. From our 
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FIG. 5. (a): Open symbols indicate the numerical results for 
the local heat flux averaged over time and angular position 
as function of Ra at different radial positions r/L. The solid 
lines show the experimental data of Shang et al. [T8|. (b): 
The scaling exponent for the time- and angle-averaged heat 
flux as function of the radial position r/L for the simulations, 
experiment [18], and theory |11| . 



numerical calculations of the kinetic side wall BL thick- 
ness, as identified by the location of the maximal velocity 
fiuctuations, we get even slightly smaller values. Obvi- 
ously, there is no thermal BL at the side walls due to the 
adiabatic boundary conditions. 

The first striking observation from Figure [2] is that 
the absolute value of the heat flux is much larger at the 
side walls (solid lines) as compared to the center (dashed 
lines). This is consistent with the expectation that most 
of the heat is transported by the large-scale convection 
roll and of course it is well known [18 . In addition, we 
observe that the PDF close to the side walls has a positive 
skewness due to rising and falling of the warm and cold 
plumes. We observe some events with a heat flux as high 
as fifteen times the average. The PDF of the local heat 
fiux in the cell center also has a positive skewness, which 
indicates that plumes can travel through this region as 
well. The inset of Figure [2] shows j as a function of the di- 
mensionless time t close to the cell center (r/L=0.06) and 
to the side wall (r/L=0.45) for Ra — 10® and again re- 
fiects the presence of much stronger heat transport events 
near the side wall than in the center. 

One of the main features of turbulence is the small 
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FIG. 6. Panel (a) shows the Reynolds number based on the 
rms velocity, Re! , and panel (b) the normalized rms temper- 
ature Q' as function of Ra\ both for different radial positions. 

scale intermittency that is measured as departure from a 
Gaussian character of the PDF, mainly the tails and the 
peakedness. This can be quantified by calculating the 
flatness of the PDF. For strongly intermittent signals 
however the integral of j'* x PDF{j) defining the flatness 
may not converge. To examine this issue we calculate the 
angular average of this quantity at r/L = 0.06 and at 
r/L = 0.45 on the mid-plane. While at the side wall this 
quantity decays for large | j \ sufficiently fast, see Fig- 
ure [3] (this feature is more evident when data is plotted 
on a linear rather than log-scale), and thus permits the 
calculation of the flatness (showing strong intermittency, 
F4 « 11, inset of Figure [s]), in the center the intermit- 
tency is so strong that no convergence for the flatness 
can be achieved. 

By rescaling the heat flux j with its standard devia- 
tion jrrns, the zcro-mcan PDF for the normalized heat 
flux j* = (j — j) firms shows Universality near the side 
wall (see Figure [4^). The tails for the rescaled PDFs 
are shorter at the side walls compared to those at the 
center. This indicates relatively fewer plumes carrying a 
large heat flux at the side walls in contrast to relatively 
more plumes carrying a smaller heat flux at the center, 
and again underlines the extremely strong intermittency 
of the heat flux in the center (see Figure ^p) . 

Figure [5^ shows the time- and angularly- averaged heat 
flux as a function of Ra at different radial positions on 
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FIG. 7. The variations in the velocity and temperature scal- 
ing exponents with Ra are shown as functions of the radial 
position. 

the mid-plane. Figure [5]d shows that the corresponding 
scaling exponent jj as a function of the radial position 
decreases monotonically from 0.43 near the axis to 0.29 
close to the side wall. The measurements of Shang et 
al. [18j and the theoretical analysis of Grossmann and 
Lohse [TT] only made statements on the values close to 
the side wall and to the axis; these are well confirmed 
by the present results. In the Ra range considered here 
the heat transport near the side wall is an order of mag- 
nitude larger than in the center. An extrapolation of 
the power law fits j = 0.0025Ra°-^^'^° °^ obtained near 
the center and j = 0.3236i?a°'^^*° °^ valid near the wall, 
shows that these become equal for Ra « 10^^. This value 
is consistent with the prediction of the unifying theory 
[5]. Recent experiments suggests the occurrence of 
this feature in the range 10^'^ < Ra < 5 x lO^"*, whereas 
Shang et al. [TS] suggested that it happens at Ra sa lO^"*, 
based on their extrapolation. 



B. Local fluctuations 

In this section, we determine the scaling with Ra of 
the velocity and temperature fiuctuations with respect to 
their global mean at different radial locations. Since the 
global mean of the vertical velocity vanishes, to charac- 
terize the velocity fluctuations Urms it is sufficient to take 
the root mean square of the fluid velocity at each posi- 
tion and average it in the azimuthal direction. Figure [6^ 
shows Re' = ^ {Ra/ Pr)urms (equivalent to UdimL/v 
with Udun the dimensional rms velocity) as a function 
of Ra on the mid-plane. 

For the normalized temperature fluctuations we take 
the root mean square of 9' = 9 — 1/2. Figure [6}d shows 
the results for this quantity vs. Ra at different radial po- 
sitions on the mid-plane. At all radial positions, fluctu- 
ations of both velocity and temperature exhibit a power 
law dependence on Ra proportional to i?a''" and Ra'^c' , 
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TABLE II. Summary of the velocity and temperature scaling exponents 7u and 751 reported in several experimental (E) and 
theoretical (T) studies. The experiments mentioned below have been carried out in cylindrical cells, unless stated otherwise. 





Ra 


Pr 


r 


7u (center) 


79 (center) 


7ti (side wall) 


79 (side wall) 


Castaing et al. '29' (E) 


4 X 10' - 6 X 10^^ 


0.65 - 1.5 


1 


0.491 ± 0.002 


-0.147 ±0.005 






Castaing et al. 29 (T) 


4 X 10 '' - 6 X 10^^ 


0.65 - 1.5 


1 


3/7 


-1/7 






Sano ct al. 30 (E) 


10« - io"J 


0.64 - 1.4 


1 






0.485 ± 0.005 




Takcshita ct al. 31 (E) 


10" - 10« 


0.024 


1 






0.46 ± 0.02 




Ashkenazi et al. 32 , 33 (E) 


10^^ - 5 X 10^* 


27 - 190 


1 (square) 






0.43 ± 0.02 




Chavanne et al. 34 (E) 


10' - 6 X 10^^ 


0.7-4 


0.5 






0.49 




Daya and Ecke 35 (E) 


2 X 10** - 4 X 10'-^ 


5.46 


0.79 


0.5 ±0.03 


-0.10 ±0.02 






Nicmcla ct al. 36 (E) 


15 X 10" - 10^^ 


0.7 


1 






0.5 




Qiu ct al. 37 38 (E) 


10« - io"J 


5.4- 5.5 


1 


0.55 




0.46 




Lam et al. 39 (E) 


10" - 10" 


6 - 1027 


0.5-4.4 






0.495 (bottom) 




Grossman and Lohse 1111 (T) 


10" - 10^* 


~ 0.1 - 10 


1 


0.34 


-0.11 to -0.16 




-0.09 to -0.11 


Shang ct al. 18"' (E) 


10« - 10^" 


4.4 


1 


0.49 ± 0.03 


-0.14 ±0.03 


0.46 ± 0.03 


-0.24 ±0.03 


Present work 


2 X 10" - 2 X 10^ 


5.2 


1 


0.44 ± 0.01 


-0.18 ±0.01 


0.49 ±0.01 


-0.20 ±0.01 




dynamical frame 



FIG. 8. Sketch of the LSC in a cylindrical RB cell. The warm 
plumes go upwards on the right side, defined as a = 0, and 
go downwards on the left side of the cell, defined as a = vr. 
Due to the shape of the LSC the warm uprising fluid (cold 
down flowing fluid) is close to the side wall at z/ L — 0.25 and 
z/L = 0.50 {z/L = 0.75 and z/L = 0.50). 

respectively. In this Ra number regime the thermal fluc- 
tuations close to the side wall are an order of magni- 
tude larger than at the cell center due to the plumes 
that travel along the wall. Figure [7] shows that the cor- 
responding velocity scaling exponents increase smoothly 
from 7„ = 0.44 in the cell center to ju = 0.49 near the 
side wall. Figure [7] shows that the corresponding tem- 
perature scaling exponent decreases from jgi = —0.18 in 
the cell center to 79/ = —0.20 near the side wall. 

Table [TT| summarizes the data for the scaling exponents 
available in the literature and compares them with the 
present ones. There are some differences among the val- 
ues reported. This is due in part to the spatial depen- 
dence of this quantity, shown in Figure [7j but also to the 
use of different experimental techniques which measure 
somewhat different quantities. Overall, there is a general 
consistency among the data shown. The origin of the 



residual differences cannot be ascertained on the basis of 
the presently available knowledge and must await further 
research. 



IV. STATISTICS WITH RESPECT TO THE 
ORIENTATION OF THE LARGE SCALE 
CIRCULATION 

For r = 1, the flow in the cell is characterized by a 
large scale circulation (LSC) [30l l40ti42] as sketched in 
Figure |8] Most of the plumes travel in the LSC plane 
close to the side wall. In experiments the LSC orienta- 
tion can be detected with thermistors embedded in the 
side wall |43) which measure relatively higher and lower 
temperatures in the region of upflow and downflow |42j . 
Here we want to determine how the local heat flux de- 
pends on the location of measurement with respect to the 
LSC orientation plane. 

A. Determination of LSC orientation 

We know from Ref. [H] that, in the present range of 
Ra numbers, the period of the LSC is approximately 10 
dimensionless time units. On this basis, according to the 
last column of Table 1, our simulations are long enough 
to cover more than 400 LSC periods for the smallest Ra 
and 175 for the highest one. 

To determine the LSC orientation we use cosine fits of 
vertical velocity, temperature and heat flux over horizon- 
tal planes in the angular direction near the side wall at 
r/L = 0.45. For this purpose we use the information from 
the numerical probes placed uniformly in this direction 
(see section 2) [42l|45]. In performing the fits we need to 
take into account that the orientation of the LSC plane 
changes with time. To this end we have pre-processed the 
instantaneous data by filtering them by means of short- 
time moving averages of four different durations, namely 
4, 10, 20 and 50 dimensionless time units. We have fit- 
ted the results of each one of these time averages by an 
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expression of the form 

fi = fni + Af cos(m0j + cj)LSc), 



(4) 



where fi is the (short-time moving averaged) information 
provided by the i*^ probe at the angular position (f)i, /,„ 
is the angular mean value, ^/ is the amplitude of the fit 
and (l)LSC is the phase shift with respect to the reference 
frame of the computation; the integer m equals 1 or 2 
depending on the angular periodicity of the quantity / 
as will be clear below. 

An example of the difference between instantaneous 
and short-time-averaged data for the heat fiux j is shown 
in the top panel of Figure [o] for Ra = 2 x 10^ on the mid- 
plane of the cylinder. The dashed line is the result of 
the fit to a portion of the instantaneous data averaged 
for 20 dimensionless time units. The bottom panel of 
Figure [9] shows the filtered data and their cosine fit for 
the vertical velocity for the same positions and interval of 
time used for the top panel. The maximum of the curve 
is marked with a red cross and identifies the position of 
the LSC plane. A graph showing (j)LSC vs. time as cal- 
culated using the four different averaging times is shown 
in Figure [TO] As expected, the fiuctuations in the posi- 
tion of the LSC plane are somewhat greater when short 
averaging times are used. 

Once the position of the LSC plane is known, all the 
short time averages can be referred to it by using a new 
angular co-ordinate a = (f) — (pLSC- With this construc- 
tion, the same value of a in each filtered segment of data 
identifies the same position relative to the LSC and jus- 
tifies taking the average over all the filtered data sets. 
Such averages for the vertical velocity and heat flux at 
r/L = 0.45 on the mid-plane are shown in Figure 11 A 



comparison of the two panels in this figure shows that 
the heat flux has double the periodicity of the velocity 
(i.e. the integer m in equation Q equals 2 as opposed 
to m = 1 for the velocity), because the heat flux is en- 
hanced in correspondence of both the upward and down- 
ward moving streams of the LSC. The local heat flux is 
lower in the areas where the LSC is not present. 

The panel with the heat flux in Figure [TT1 actually con- 
tains four different sets of symbols, each one denoting the 
results obtained with a different averaging time. The vir- 
tual coincidence of these four sets demonstrates the ro- 
bustness of the results with respect to the filtering time 
used. 



In Figure 12 the amplitudes of the cosine fits for tem- 
perature, Ag, axial velocity, A^^, and normalized heat 
fiux, Aj are shown as functions of Ra at mid-height near 
the side wall, r/L = 0.45. Interestingly, these ampli- 
tudes have a power-law dependency on Ra, with scaling 
exponents -0.250 ± 0.010, 0.020 ± 0.005, 0.019 ± 0.005 
for temperature, axial velocity and local heat flux respec- 
tively. 

Figure [l3^ shows the time-averaged local heat flux as a 
function of Ra at different a at mid-height near the side 
wall. From the flgure it is clear that the local heat flux 
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FIG. 9. Angular dependencies of heat flux and vertical ve- 
locity on the mid-plane at r/L = 0.45 near the side wall for 
Ra = 2 X 10*: (a) instantaneous local heat flux (green points), 
(b) vertical velocity (blue points). The orientation of the LSC 
plane is determined by using the flltered time average of 20 
dimensionless time units on vertical velocity data and the best 
cosine flt is shown in black solid line. The red cross indicates 
the maximum of the cosine flt, which we take as the LSC 
orientation. 
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FIG. 10. LSC orientation as function of time for Ra — 2x 10*. 
The time averaging that is applied to the signal before the 
analysis is 1 (blue-solid), 4 (red-dashdot), 20 (black-solid) and 
50 (magenta-dash) dimensionless time units. 



has a power law dependence on Ra. In agreement with 



the data in Figure 12 we find that the local heat fiux at 
the LSC orientation increases faster than in the regions 
away from LSC. This is also revealed when the local heat 
fiux scaling exponent as function of a, see Figure [Tsj j, is 
considered. 



Figure 14 shows the variation of the time-averaged 
local heat flux with a at three different heights. For 
z/L — 0.25 larger values of j occur near a = 0, and 
lower values near a = tt. A similar picture shifted by tt 
is found for z/L = 0.75, with higher values near a = -k 
and lower values near a = 0. At z/L = 0.50, on the 
other hand, the levels at a = and a = tt are compara- 
ble. These results suggest that the plane of the LSC is 
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0.5 0.75 

a /{2tt ) 

FIG. 11. Time averaged quantities relative to the orientation 
of the LSC plane near the side walls {r/L = 0.45) at mid- 
height: (a) vertical velocity averaged over 20 dimensionless 
time units; (b) local heat flux j{a) averaged over 4 (circle- 
blue), 10 (square-red), 20 (diamond-black) and 50 (triangle- 
magenta) dimensionless time units. The best cosine fit, with 
a period of Air which corresponds to m = 2 in Q , is shown by 
dashes. The mean heat flux (thin-dash dot) and the ampli- 
tude of variation Aj are also shown. The hot plumes originate 
is in the neighborhood of a = 0. 
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FIG. 12. The amplitude variation of cosine fits near the 
side walls, r/L = 0.45, at the mid-height as function of Ra. 
The data are shown for temperature (circle-red), axial veloc- 
ity (square-black) and normalized heat flux (diamond-blue). 
Here | j{a) \m,ean is the arithmetic mean heat flux in the az- 
imuthal direction for a given Ra. The dashed lines indicates 
the power law flts for the data. 

tilted as sketched in Figure [8| This feature has also been 
found in experiments of Funfschilling et al. |46j . 



V. SUMMARY AND CONCLUSIONS 

To summarize, we investigated numerically the scaling 
of the local heat flux in Rayleigh-Benard convection of a 
fluid vi^ith Pr — 5.2 (appropriate for water at 32 °C) for 
2 X 10^ < Ra < 2 X 10^ in an unit aspect ratio cylinder. 
In this Ra number regime the local heat flux is larger 
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FIG. 13. (a): Scaling of the time averaged local heat flux with 
Ra near the side walls, r/L — 0.45, at mid-height. The sym- 
bols indicate measurements taken at different a. The dashed 
lines are the power law flts to the data, (b) : Scaling exponent 
for the heat flux 7j (a) relative to the LSC plane near the side 
wall, r/L = 0.45 at z/L = 0.25 (triangles-red), z/L ^ 0.50 
(square- black), and z/L — 0.75 (diamond- blue). The straight 
lines indicate the arithmetic mean values of the scaling expo- 
nents. 



close to the side wall than on the axis. The local heat 
flux Uz{T — To) is a positive quantity both when fluid 
warmer than the average temperature Tq rises and fluid 
colder than Tq sinks. The PDFs of the local heat flux 
have a positive skewness due to the dominance of plume 
transport in this Ra range. On the mid-plane, the scaling 
exponents of the local heat flux with Ra near the axis and 
close to the side wall agree well with the measurements 
of Shang et al. [TS] and the predictions of Grossmann 
and Lohse [11]. Here we have shown that these scaling 
exponents decrease monotonically with position r from 
0.43 near the axis to 0.29 close to the side wall. The 
scaling exponent for the Reynolds number based on the 
rms velocity depends on the radial position as well, with a 
value of 0.44 near the axis and 0.49 close to the side wall. 
For the scaling exponents of the temperature fluctuations 
we find —0.18 and —0.20 respectively. We showed the 
marked effect of the LSC which causes local heat fluxes 
more than twice as large in its plane than at 90° from it. 
This effect becomes stronger for high Ra. 
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